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ABSTRACT 

Hybrid Vlasov-Maxwell simulations are employed to investigate the role of kinetic effects in a two- 
dimensional turbulent multi-ion plasma, composed of protons, alpha particles and fluid electrons. In 
the typical conditions of the solar-wind environment, and in situations of decaying turbulence, the 
numerical results show that the velocity distribution functions of both ion species depart from the typ- 
ical configuration of thermal equilibrium. These non-Maxwellian features are quantified through the 
statistical analysis of the temperature anisotropy, for both protons and alpha particles, in the reference 
frame given by the local magnetic field. Anisotropy is found to be higher in regions of high magnetic 
stress. Both ion species manifest a preferentially perpendicular heating, although the anisotropy is 
more pronounced for the alpha particles, according with solar wind observations. Anisotropy of the 
alpha particle, moreover, is correlated to the proton anisotropy, and also depends on the local dif- 
ferential flow between the two species. Evident distortions of the particle distribution functions are 
present, with the production of bumps along the direction of the local magnetic field. The physical 
phenomenology recovered in these numerical simulations reproduces very common measurements in 
the turbulent solar wind, suggesting that the multi-ion Vlasov model constitutes a valid approach to 
the understanding of the nature of complex kinetic effects in astrophysical plasmas. 



1. INTRODUCTION 

Astrophysical plasmas are generally in a fully turbulent 
regime. In particular, the solar wind can be considered 
as a natural laboratory for studying physical processes 
of plasma turbulence, whose dynamic scales cannot be 
achieved in laboratory experiments. Collisions in the 
solar wind are very rare: the particle mean free path 
is comparable to (or larger than) the system size (~ 1 
AU). These properties suggest that small-scale mecha- 
nisms should be more complex than one would expect 
I from a fluid (viscous) description. 

Spacecraft measurements generally reveal that the 
electromagnetic field fiuctuations in t he solar wind are in 
a stat e of fully -developed turbulence (jBruno &: Carbon^ 
! 120051: lMarschl l200l . The power spectrum of the fluc- 
tuating fields manifests a behavior reminiscent of the 
■ power-law for fluid turbulence (jKolmogorovll 1 9411 ) . 

where k here is the wavenumber obtained applying the 
Taylor hypothesis. This inertial range turbulence ex- 
] tends to smaller spatial scales down to a range of wave- 
lengths where kinetic effects dominate the plasma dy- 
namics. At these scales, different physical processes 
come into play, leading to spectral shape changes. The 
first clear spectral change appears at scales such as 
the ion inertial length and/or the ion Larmor radius 
(iLeamon et al.l [20001 : iBourouaine et alll20l3 iBale et all 
l2005f ). where the spectr um of the magnetic field gets 
steeper (jBale et al.ll2005[ ). At the range of spatial lengths 
of the order of the electron kinetic scales the interpre- 
tation of the solar-wind observations is still controver- 
sial. Two different scenarios have been recently pictured: 
a seco nd spectral break wi th an additional power-law 
range (ISahra oui et al.ll2010D and an exponential cut-off 
(jAlexa ndrova et al.ll2009[ )~tliat instead marks the end of 
self-similarity. In both theoretical support from 



self-consistent, fully nonlinear, Vlasov models is needed 
for the interpretation of this complex phenomenology. 

This general picture of astrophysical turbulence be- 
comes more complicated because of the multi-component 
nature of the solar wind. The interplanetary medium, 
although predominantly constituted of protons, is also 
made of a finite amount of doubled ionized helium (al- 
pha particles), together with a few percentage of hea v- 
ier ions. Several observations (M arsch et al.l Il982ai rbl: 
'Kasp er et al.l [2008^ have shown that heavier ions are 
heated and accelerated preferentially as compared to pro- 
tons and electrons. Moreover, in recent analysis per- 
formed on solar- wind data from the Helios spacecraft, 
the link between the signatures of kinetic effects and 
some important parameters of heavy ions, such as rel- 
ative speed, temperature ratio a nd anisotro py, has been 
investigated (jBourouaine et al.ll2010ll2011al fa). In these 
works the authors pointed out that more significant 
anisotropies and non-Maxwellian features are detected 
for alpha particles distribution functions with respect to 
protons. The evolution of the velocity distribution func- 
tions in the solar wind, and the production of kinetic 
signatures such as heating and temperature anisotropies 
represen t nowadays some of the key issues of plasma 
physics (jOsman et al.ll201lL [20T1) . 

The two points discussed above, that the kinetic scales 
have a determinant effect in shaping the turbulent spec- 
tra and that the role of secondary ions cannot be ne- 
glected, suggest that a multi-scale and multi-species self- 
consistent Vlasov treatment of the turbulent solar wind 
is required. The detection of the main kinetic processes 
that, acting on small scales, are responsible for the en- 
ergy dissipation and the heating production, represents 
today one of the main challenges for the plasma scientific 
community. A kinetic description of coUisionless plasma 
turbulence offers the powerful opportunity of giving im- 
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portant insights for the interpretation of "in situ" satel- 
lite measurements in the solar wind. In this context, 
an indispensable and crucial support to investigate the 
complexity of solar- wind physics is represented by kinetic 
numerical simulations. 

The most widely adopted numerical description of the 
kinetic plasma physics is represented by the Particle in 
Cell (PIC) methods. The Lagrangian PIC approach is 
based on the integration of the equations of motion of 
a large number of macro-particles under the effects of 
self-consistent electromagnetic fields. The PIC simu- 
lations have been extensively used for the description 
of the kinetic dynamics of space plas mas, addressing 
the p roblem of wave-pa rticle interaction (lAraneda et al.l 
I2008D . particle heating fAra neda et a l."2009^ and turbu 

2008; Parashar ct al 
Parashar et al.ll201ll 



lence |Garv ct al. 2008; Sait o et al. 
2010"; Camoorcalc & BurEcss|l20lIF 
Markovskii &: Vasque z 2002|). Nowadays, thanks to the 
fast technological development of the computational re- 
sources, the Eulerian approach for the numerical solution 
of the Vlasov equation has b ecome accessible as aii alter- 
native to the PIC approach CMange nev et al.ll2003 ). The 
Eulerian Vlasov simulations, in which the time evolution 
of the particle distribution function is followed numer- 
ically in a discretized phase space domain, are signifi- 
cantly more demanding from the computational point of 
view compared to the PIC simulations. Nevertheless, at 
variance with the PIC methods, the Vlasov algorithms 
are mostly noise-free and do not introduce additional 
(unphysical) heating due to the particle noise. This point 
can be crucial when dealing with the numerical descrip- 
tion of the short wavelengths part of the turbulent cas- 
cade, where the energy level of the fluctuations is typi- 
cally very low and the statistical noise introduced by the 
PIC calculations can possibly cover the physical infor- 
mation. 

Re cently ()Mangenev et al.ll2002t [Valentini et al.ll2005l 
120071) . an Eulerian hybrid Vlasov- Maxwell code (HVM 
hereafter) has been developed. This algorithm integrates 
numerically the Vlasov equation in phase space coupled 
to the Maxwell equations for the electromagnetic fields. 
The Vlasov equation is solved for the proton species, 
while electrons are treated as a fluid (jValentini et al.l 
|2007[) . The HMV code has been extensively used for 
the analysis of the kinetic effec ts during the evolu- 
tion of the solar- wind cascade fValentini ct al. 2008'; 
iValentini fc Veltril 120091: IValentin i ct al. 2010, 2011). 
Furthermore, IServidio et al.l(|2012D made use of the HMV 
algorithm to investigate the role of local kinetic effects in 
plasma turbulence in a 2D-3V phase space configuration 
(two dimensions in physical space, and three in veloc- 
ity space). It has been shown that, nearby the region of 
strong magnetic activity, the proton distribution function 
is deformed by kinetic effects displaying significant non- 
Maxwellian features. Moreover, in these regions charac- 
terized by high magnetic stress, reconnection events can 
occur locally, as it has been also describe d in both PIC 
and Magnetohydrodynamic simulations (jServidio et al.l 
[200a Prake ct'al. 2010). 

As discussed previously, in the solar wind description 
it is important to take account of its multi-c omponent 
nature. For this purpose. iPerrone et al.l ()201lD have pro- 
posed an update version of the HVM code, which in- 
cludes the kinetic dynamics of heavy ions. In particular. 
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Fig. 1. — Time evolution of the average out-of-plane squared 
current density (jf). 

the authors examined the effects produced by the pres- 
ence of alpha particles in the evolution of the solar-wind 
turbulent cascade in the direction parallel to the ambient 
magnetic field in a 1D-3V phase space configuration. 

In the present paper we discuss the results obtained 
by the "multi-component" version of the HVM code in 
a 2D-3V phase space configuration. In order to study 
the turbulent activity in the presence of alpha parti- 
cles, we adopt an ap proach similar to the one used by 
IServidio et al.l (|2012| ) , extending the above results to the 
more realistic multi-ion treatment. As we will discuss in 
the following, our numerical simulations reproduce sev- 
eral features commonly observed in space plasmas. As 
the result of the turbulent cascade, coherent structures 
appear in the system, with regions of high magnetic 
stress where reconnection locally occurs; in correspon- 
dence of these regions the particle distribution functions 
are heavily distorted, exhibiting a significant tempera- 
ture anisotropy. This effect is more evident for alpha 
particles than for proton s, consistent with recent obser - 
vations in the solar wind (jBourouaine et aO2010l[2011b( ). 

2. NUMERICAL MODEL 

We simulate a collisionless and magnetized multi- 
species (electrons, protons and alpha particles) turbulent 
plasma through the use of the HVM code. Within this 
HVM model, the Vlasov equation for pro ton (/„) and al- 
pha particle (/„.) distribu tion functions (jValentini et al.l 
l2007t [Perrone et al.ll201ll ) is integrated numerically in a 
2D-3V phase space domain. Electrons are treated as a 
fluid and a generalized Ohm equation, where a resistive 
term has been added as a standard numerical Laplacian 
dissipation, is considered. The dimensionless HVM equa- 
tions are given by 



f +v.^ + (E + vxB).fi^O, 
at or av 



E 

'dt 



-(Ue X B) - — VPe 



= -V X E. 



(1) 

(2) 
(3) 



The quasi-neutrality approximation = ZiUi is con- 
sidered, where rie and are the electron and ion den- 
sities respectively and Zi the ion charge number (the 
subscript i = p,a stands for protons and alpha particles 
respectively). In Eq. ([T]), /^(r, v, t) is the ion distribution 
function, E(r, t) and B(r, t) are the electric and magnetic 
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fields. In Eq. the electron bulk velocity is defined 
as Ue = {J2i ZiUiUi — y X B)/ne, where the ion bulk ve- 
locities Ui are evaluated as first order velocity moments 
of the ion distribution function; finally, j = V x B repre- 
sents the total current density. An isothermal equation 
of state for the electron pressure closes the system. 

In Eqs. ([I])-® time is scaled by the inverse proton- 
cyclotron frequency il~p, velocities by the Alfven speed 
Va, lengths by the proton skin depth dp = Va/^cp and 
masses by the proton mass nip. From now on, all the 
physical quantities will be expressed in units of the char- 
acteristic parameters listed above. A small value for 
the resistivity (77 = 2 x 10~^) has been chosen in or- 
der to achieve relatively high Reynolds numbers and to 
remove any spurious numerical effects due to the pres- 
ence of strong current sheets. The initial equilibrium 
consists of a plasma composed of kinetic protons and al- 
pha particles, with Maxwellian velocity distributions and 
homogeneous densities (no,p and rio^o- respectively), and 
fluid electrons embedded in a background magnetic field 
Bo = BqBz- The plasma dynamics and the development 
of turbulence are investigated in a double periodic (x, y) 
domain perpendicular to Bq, where the total magnetic 
field can be written as B = Bq + B^. 

The equilibrium configuration is perturbed by a 2D 
spectrum of fluctuations for the magnetic and proton ve- 
locity fields (alpha particles have zero initial bulk veloc- 
ity). We inject energy with random phases and wave 
numbers in the range 0.066 < fc < 2, where k = 27rm/L, 
with 2 < TO < 6 and L being the box size in each spatial 
direction. The rms of the initial magnetic perturbations 
is 6B/B0 ~ 0.3. Neither density disturbances nor paral- 
lel variances are imposed at t = 0. The proton plasma 
beta is (3p = 2vth,p'^/V2 = 2 (where Vth,p = y/Tp/nip = 1 
is the proton thermal speed) and the electron to proton 
temperature ratio is Te/Tp — 1. For the alpha parti- 
cles we set Za = 2; to^/top = 4, ?T.o,a/"-o,p = 5% and 
Ta/Tp = 1. With this choice, the alpha particle thermal 
speed is vth,a = vth,p/'2- 

The system size in the spatial domain is i = 27r x 2Qdp 
in both X and y directions, while the limits of the velocity 
domain for both ion species are fixed at Vmax,i — i5vth,i 
in each velocity direction. In these simulations, we use 
512^ grid-points in the two-dimensional spatial domain 
and 61^ and 31'^ grid-points in proton and alpha par- 
ticle three-dimensional velocity domains, respectively. 
We point out that in the Ohm equation for the elec- 
tric field we have neglected the electron inertia terms. 
These terms are in fact proportional to the squared 
electron skin depth (which in scaled units is given by 
— TOe / TOp) , then cannot be adequately resolved within 
the discretization of our simulations. The time step 
At has been chosen in such a way that the Courant- 
Friedrichs-Lewy condition for th e numerical s tability of 
the Vlasov algorithm is satisfied (jPevret fc Tavlorn 19861 ). 

3. NUMERICAL RESULTS 

We numerically study the kinetic evolution of protons 
and alpha particles in a situation of decaying turbulence. 
We expect that kinetic effects develop simultaneously 
with magnetic fluctuations and shears, the latter play- 
ing a fundamental role in the production of interesting 
features such as particle acceleration, heating, tempera- 
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Fig. 2. — Contour plots (shaded colors) of the out-of plane current 

densities: jz (a), ji'^"^ (b), ji^'^ (c), and j'^"^ (d). The isolines 
of the magnetic potential Az are indicated by black/white lines. 
In panel (a), the positions of the X-points (regions of magnetic 
reconnection) are indicated by black thick crosses. 
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Fig. 3. — In the top panel we show the PDFs of the different 

(e) (v) 

current densities: jl red sohd-cross hne, blue dashed-square 

(a) 

line, ji green dashed-circle line, and jz black dashed-triangle 
line. In the bottom panel, the PDFs of the standardized variables, 
obtained by subtracting the average and normalizing to the respec- 
tive rms value, are reported. The light-blue dot-dashed line is the 
Gaussian fit. 



ture anisotropy, wave-particle like interactions, and gen- 
eration of beams in the ion distribution function. 

As in the fluid counterpart, large scale fluctuations pro- 
duce a turbulent cascade toward small scales (high fc's). 
In analogy with fluid n iodels (MHD, Hall MHD etc.) 
of decaying turbulence (jMininni fc PouQuetl [20091) ■ it is 
possible to identify an instant of time at which the tur- 
bulent activity reaches its maximum value. Since the 
current density is proportional to the level of small-scale 
gradients, a good indicator of the level of turbulent ac- 
tivity is represented by the average out-of-plane squared 
current density (j^), whose time evolution is shown in 
Figure [H Evidently, a.t t — t* AO {jl) attains its 
maximum value. This is the characteristic time at which 
decaying turbulence shares many statistical similarities 
with steady state (driven) turbulence, and, at this time, 
we perform our analysis. 

The turbulent activity leads to the generation of co- 
herent structures. Vortices and current sheets appear 
in the contour plots of Figure [51 where shaded contours 
of the out-of-plane total current density jz [panel (a)] 
are represented. In the same figure, for the descending 

panels, the separate contributions ji of each species, 
namely j| (electron current), j^ (proton current), and 
(alpha-particle current) are also shown. The contour 
lines in each plot represent the magnetic potential Az 
of the inplane magnetic field (B_l = VAz x e^). The 
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Fig. 4. — Power spectra of rip (green-square line), Up (purple- 
triangle line), B (black-solid line) and E (red-dashed line). 

different colors (black/white) of the Az contour lines in- 
dicate different directions of rotation of the vortices. At 
a careful analysis of the four plots of Figure [H a certain 
correlation seems to exist between the proton and the 
alpha-particle current densities, revealing that the local 
small-scale structures of the two ion species behave in a 
similar way. Moreover, the electron flow generates struc- 
tures at smaller scales, manifesting a more intermittent 
and bursty character. The coherent structures visible in 
the four plots of Figure [D are not static, but evolve in 
time interacting nonlinearly among each others. More- 
over, in between the islands, the current becomes very 
intense and magnetic reconnection events locally occur 
at the X points of A^, indicated in the contour plot of 
panel (a) by black crosses. The presence of these high 
magnetic stress regions is a signature of the intermittent 
nature of the magnetic field, and this affects also the 
patchiness of the parallel and perpendicular heating (see 
below). 

In the top panel of Figure [3| we report the probabil- 
ity distribution functions (PDF) of j'f^ (red solid-cross 

line), j'f^ (blue dashed-square line), j'f^ (green dashed- 
circle line) , and jz (black dashed-triangle line) . This plot 
clearly indicates that there is a certain ordering in the 
maximum values of the achieved current. The main con- 
tribution to the total current seems to come from the 
electrons and the protons, that develop the most intense 
bursty events. In contrast, the alpha-particle current 
structures are smoother and are concentrated on larger 
scales, as can be noticed from Figure [2] In the same fig- 
ure (bottom panel) , we report the PDF of the standard- 
ized variables obtained by subtracting the average and 
normalizing to the respective rms value. The Gaussian 
fit is also plotted (light-blue dot-dashed line) as refer- 

(e) 

ence. The currents jz and jz are highly non-Gaussian 
distributed, because they are related to the increments 
(gradients) of the magnetic field (and electron flows are 
essentially frozen-in). The proton and alpha particle con- 
tributions, on the other hand, behave more like Gaussian 
variables, since they are related to primitive variables of 
turbulence such as velocities and densities, and they do 
not capture high order statistics. 
In order to quantify the turbulence, the power spec- 
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Fig. 5. — Panel (a): power spectra of the proton bulk 

velocity (i = p, black-solid line) and of the alpha-particle bulk 
velocity (i = a, red-dashed line); panel (b): power spectra l-^V^^'p 
of the proton density (i = p, black-solid line) and of the alpha- 
particle density (i = a, red-dashed line) 



tra of density and bulk velocity for both protons and 
alpha particles and of magnetic and electric fields have 
been computed. In Figure 2] we show the power spec- 
tra of Up (green-square line), Up (purple-triangle line), 
B (black-solid line) and E (red-dashed line). The Kol- 
mogorov expectation k~^^^ (black-dashed line) has been 
plotted in Figure U as a reference. These omnidirectional 
power spectra reveal several interesting features, many 
of them recovered also in solar-wind spacecraft observa- 
tions. In fact, the large scale activity is incompressible 
and the Alfvenic correlation between magnetic and ve- 
locity fluctuations is broken at the proton skin depth 
(vertical black dashed line). Moreover, the electric ac- 
tivity (red-dashed line) at higher wavenumbers is sig- 
nificantly more intense than the magnetic one (bl ack- 
solidline) (jBruno fc Carbonel[2005l: iBale et al.ll2"005[) . It 
is worth to point out that the power spectra displayed in 
Figure |4] present no significant differences with respect 
to the same spectra obtained through HV M simulations 
without alpha particles [see Figure 1(b) in lServidio et alJ 
(|2012D ]. meaning that the presence of alpha particles does 
not significantly affect the dynamical evolution of the 
turbulent cascade. To make a direct comparison of the 
dynamical evolution of the two ion species, in Figure [5] 
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Fig. 6. — Top panel: PDF of the temperature anisotropy Ai = 

T^*'/T|J'' of protons (a) and alpha-particles (b) at four different 

times in the system evolution: t = 1 black line, t = 21 purple 
line, i = 34 red line, f = 40 blue line. Bottom panel: comparison 
between the proton (Ap, black line) and the alpha- particle {Aa, red 
line) anisotropics is reported at the peak of the turbulent activity, 
i.e., a.t t = 40. 



we show the velocity spectra 
normalized density spectra 



W|2 



7V*;'^|2 [panel (b)] for pro- 



[panel (a)] and the 



tons [i = p, black-solid lines) and for alpha particle 
{i = a, red-dashed lines). The density spectra for pro- 
tons and alpha particles are normalized to no,p and no,Q 
respectively. While the velocity spectra [panel (a)] of the 
two ion species do not display significantly different fea- 
tures, we notice that the alpha particles contribution to 
the density spectra [panel (b)] is lower than the proton 
one for wavenumbers higher than the proton skin depth 
wavenumber. This behavior is possibly related to the 
fact that the alpha particles are heavier than protons, so 
their inertia does not allow to follow the field fluctuations 
at smaller scales. 

At this point it is important to investigate the link 
between the turbulent behavior observed in the plasma 
and the generation of non-Maxwellian features in velocity 
space. For this purpose, the ion temperature anisotropy 

for each species A, = t|*VT||'\ defined as the ratio be- 
tween the perpendicular and the parallel temperature 
with respect to the local magnetic field, has been com- 
puted. Our initial condition has been set up in such a way 
to have spatially isotropic temperatures for both the ion 
species at t — 0. Nevertheless, during the development 
of turbulence the temperatures do not remain spatially 
isotropic but present local enhancements and depressions 
nearby the regions of high magnetic stress (not shown), 
as already found for the protons in .Servidio et al.. (2012J . 
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Fig. 7. — Joint probability distribution function of proton and 
alpha-particle temperature anisotropy. This two-dimensional PDF 
shows a correlation between the anisotropy in the two species. The 
best fit Ap ~ Aa, with f = 0.22, is also reported with a thick solid 
line. This behavior is in good agreement with solar wind data, as it 
can be seen comparing this plot with Figure l-(b) in lMaruca et al.1 

poll ). 



Figure |6] shows the probabihty distribution function 
(PDF) of the temperature anisotropy for protons Ap 
[panel (a)] and alpha particles Aa [panel (b)] at four dif- 
ferent times in the simulation. In the early stage of the 
system evolution (t = 1 black line), the PDFs are picked 
around Ap ~ Aa = 1, meaning that the simulation starts 
with an isotropic configuration. During the evolution of 
the system {t — 21, purple line; t = 34, red line; t — 40, 
blue line) the PDFs elongate in the parallel {Ai < 1) 
and in the perpendicular {Ai > 1) direction, displaying a 
strong anisotropic behavior. It is worth noting that the 
statistical behavior of the anisotropics saturates already 
at t ~ 20V,~p. Regardless of the particular ion species, 
the anisotropy preferentially manifests itself along the 
perpendicular direction, an evid ence commonly detected 
in the solar- wind observations (jBourouaine et all 120101 
l2011bD . However, alpha particles are more anisotropic 
than protons, as more evidently shown in Figure [6] (bot- 
tom plot) where we have directly compared the PDFs of 
the two ion species (Ap black line and Aa red line) when 
the peak of the nonlinear activity is reached {t = AOiljp). 

A question now naturally occurs: are these patchy 
anisotropics correlated? Any correlation between Aa and 
Ap may reveal that simultaneous kinetic instabilities lo- 
cally occur, modulated by the ambient magnetic field, or 
that an instability for a given species may influence the 
other, and vice- versa. In Figure [7] we analyze the cor- 
relation between protons and alpha particle temperature 
anisotropy, showing the joint PDF. Although most of the 
events are concentrated at Ap — Aa — 1 (isotropy), and 
are broadly scattered because of turbulence, this joint 
distribution suggests that there is a clear monotonic de- 
pendency between alpha and proton anisotropics. The 
shape of this distribution is in good agreement with so- 
lar wind data, a s it can be seen comp aring our Figure [7] 
with Figu re 1 o f ^Maruca et al.l (|2012[ ). Moreover, analo- 
gously to iMaruca et al- (,2013) . we fitted the above dis- 
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Fig. 8. — Anisotropy of the alpha-particle temperatures binned 
as a function of the differential speed Vap = jup — Ua| (in Alfvenic 
units). This behavior is in go od agreement with solar w ind obser- 
vations (cfr. with Figure 4 of IIBourouaine et ani201lEfl ). 



tri bution with An = obtaining C — 0.22 (note that 
in IMaruca et al.l (|2012l ) the authors obtained C — 0.25). 
These results suggest that the correlation between proton 
and alpha particle kinetic effects, commonly observed in 
the solar wind, may be the result of an active turbulent 
cascade, where kinetic instabilities are locally activated 
and modulated by the ambient magnetic field. 

In a multi-ion plasma another source of instability is 
represen ted by the diffe rential fl ow between different ion 
species (|Ganl[l993; Garv et al.i r2006,. ,2008,). We found 
that the temperature anisotropy for the alpha particles 
shows a certain correlation to the their drift velocity 
Vap = |up — Uq,| with respect to protons. This is visi- 
ble in Figure m where we report Aa as a function of Vap'- 
the temperature anisotropy increases with increasing rel- 
ative flow speeds (in Alfvenic units), up to Vap ~ 0.5. 

These results are aga in in good agreem e nt with some 
observational analyses. iBourouaine et al.l ()2011bD stud- 
ied correlations of temperature anisotropics and differ- 
ential ion speed in the solar-wind measurements from 
the Helios spacecraft; for the case of the alpha parti- 
cles, they found that Aa increases as the ion differential 
speed stays below about 0.5Va- Beyond this value Aa 
becomes roughly constant, until Vap exceeds a value of 
about 0.7Va, but then it decreases towards a value below 
unity when Vap — Va (not reached in our syste m). Com- 
paring our Figure |H] with Figure 4 of Bourou aine et aP 
(|2011bD . we flnd a very good correspondence. However, 
it is worth to point out that in a different data analy- 
sis of Advance d Composition Exp lorer (ACE) solar wind 
observations, ()Kasper et al.l [20081 ) found that the alpha 
temperature anisotropy is monotonically decreasing with 
increasing alpha particle to proton relative speed in the 
range < Vap ^0.5. It is also worth noting that while 
these studies are carried out on years of solar wind data, 
that detect different plasmas with different parameters, 
homogeneities, large scale effects and so on, in our case, 
these phenomena are the genuine result of a turbulent 
and statistically homogeneous cascade. 

To conclude our study, we discuss a few examples of 
the effects of turbulence on the velocity distributions of 
alpha particles. In Figure [HI we show the iso-surfaces of 
the alpha-particle velocity distribution at two distinct lo- 
cations in physical space, at which Aa > 1 [(a)-(b)], and 
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Fig. 9. — Iso-surfaccs of the alpha particle velocity distribution 
fa{r* ,Vx,Vy,Vz), at two different spatial locations: namely in re- 
gions where the distribution function shows anisotropy Aa >1 (a)- 
(b) and Aa <1 (c)-(d). The direction of the local magnetic field 
(red tube) is also reported together with the princip al axis (blue 
tube) of the velocity distribution (see lServidio et al.l lj2012l ')'). 

Aa < 1 [(c)-(d)]. In the same figure, we also report the 
direction of the local magnetic field (red tube) and the 
principal axis (blue tube) of the velocity distribution, 
evaluated from the stress tensor in the minimum va ri- 
ance frame (for more details see IServidio et all ([Ml V 
The alpha-particle velocity distribution appears strongly 
affected by turbulence and modulated by the local mag- 
netic field topology, manifesting both kinds of anisotropy; 
moreover, as it is clear from Figure [9l the principal axis of 
the velocity distribution can be both aligned or perpen- 
dicular to the local magnetic field. Another interesting 
feature that can be appreciated in the plots of Figure 
[3] is the local formation of bubbles in the velocity dis- 
tribution along the direction of the local magnetic field 
that resemble the characteristic longitudinal beams of ac- 
celera ted particles common ly observed in the solar wind 
data (Marsch et al.' '1982a''F) an d in 1D-3V HVM sim- 
ulations (Valcntini et al. 2008; Perrone et all 1201 In 
these numerical papers a possible mechanism responsi- 
ble for the generation of these beams of accelerated par- 
ticles has been identified in the resonant interaction of 
particles with a newly identified branch of electrostati c 
fluctuations called ion-bulk waves (jValentini et al.ll2011h . 
This resonant interaction resul ts more efficient for pro- 
tons than for alpha particles (jPerrone "eFalllmil). m 
agree ment with recent solar-wind data analysis (jMarschl 

4. SUMMARY AND DISCUSSION 

In this paper we discussed the numerical results of Hy- 
brid Vlasov-Maxwell simulations employed to investigate 
the role of kinetic effects in a two-dimensional turbulent 
multi-ion plasma, composed of protons, alpha particles 
and fluid electrons, in the typical conditions of the solar- 
wind environment. While we pointed out that the pres- 



ence of a small percentage of heavy ions does not affect 
the evolution of the turbulent cascade, our numerical re- 
sults clearly show that the dynamics of alpha particles at 
short spatial scales displays several interesting aspects, 
mainly consisting in the departure of the distribution 
function from the typical Maxwellian configuration. 

In a situation of decaying turbulence, coherent struc- 
tures appear, such as vortices and current sheets. In be- 
tween magnetic islands, reconnection events occur. We 
quantify the contribution to the total current due to pro- 
tons, alpha particles and electrons, finding that the most 
of the contribution comes from the electrons, while the 
current density of alpha particles is very low, due to 
their higher inertia. Moreover, we found the currents 
of protons and alpha particles are Gaussian distributed, 
while electron current densities are highly intermittent, 
since they capture the small-scale features of the mag- 
netic field. 

Power spectra for protons and for the electromagnetic 
fields are not too much affected by the presence of alpha 
particles, while at smaller scales (at k's higher than the 
inverse proton skin depth), the power spectrum of the 
alpha density has a steeper slope than that of proton 
density. 

The non-Maxwellian features generated during the 
evolution of turbulence have been quantified through the 
statistical analysis of the temperature anisotropy in the 
reference frame of the local magnetic field. The joint 
probability distribution of the temperature anisotropics 
of alpha particles and protons suggests that there is a 
clear monotonic dependency between alpha and proton 
anisotropics. This reveals the occurrence of simultaneous 
local kinetic instabilities or the infiuence of the dynami- 
cal evolution of one species on the other. In general, for 
the parameters studied here, we found that alpha parti- 
cles develop higher anisotropy than protons. 

With the purpose of analyzing the role of kinetic in- 
stabilities driven by the relative speed of the two ion 
species, we evaluated the dependence of the temperature 
anisotropy of alpha particles on their drift speed Vap with 
respect to protons. From this analysis, we found that the 
temperature anisotropy increases with increasing relative 
flow speeds up to Vap ~ 0.5 (in units of Va). Moreover, 
the velocity distributions of alpha particles create bumps 
along the local magnetic field, resem bling very com mon 
structures observed in the solar wind iMarschI (|2006i ) . 

By comparing our numerica l resu lt s to recent solar- 
wind analyses (jMaruca et al.l 120121 : iBourouaine et al.l 
l2011bD . we found a very good quantitative correspon- 
dence both for the correlation of alpha particle and pro- 
ton temperature anisotropics and for the correlation of 
alpha anisotropy and relative flow speed. However, it 
is worth noting that, while the observational studies are 
carried out on years of solar wind data, that detect plas- 
mas with very different features and in different physical 
regimes, in our case these correspondences are the gen- 
uine results of a turbulent cascade, where locally both 
the magnetic field topology and the relative motion of 
different ion fiows can be the main sources of kinetic ef- 
fects. 

The results presented in this work significantly repro- 
duce an important part of the complex phenomenology 
underlying many processes in the solar wind and sug- 
gest that a noise-free Eulerian Vlasov description of a 



8 



Per rone et al. 



multi-component collision-free plasma plays a fundamen- 
tal role in the interpretation of the observational data 
from spacecraft. 
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